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We use a correlated local Fermi gas (LFG) model, which accounts also for long 
distance corrections of the RPA type and final-state interactions, to compute the 
polarization of the final lepton in charged-current quasielastic neutrino scattering. 
The present model has been successfully used in recent studies of inclusive neutrino 
nucleus processes and muon capture. We investigate the relevance of nuclear effects 
in the particular case of r polarization in tau-neutrino induced reactions for several 
kinematics of relevance for neutrino oscillation experiments. 
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I. INTRODUCTION 

The purpose of this paper is to investigate the importance of nuclear effects on the 
final lepton polarization in neutrino-induced charged-current (CC) inclusive reactions of the 
type {I'ul). In particular we present results for the tau polarization in (^'t,t) reactions. 
The interest of these studies on neutrino-nucleus reactions lies on their implications for the 
neutrino oscillations experiments Q, i (see i and references therein for a recent review). 
Some of the experiments proposed to demonstrate the oscillation are expecting 

to detect the r production signal through the (i/^, r~) or (z7^, r"'') reactions 0,0], among 
them the CNGS project which will send a neutrino beam from CERN to the Gran Sasso 
laboratory, where the ICARUS and OPERA detectors will start taking data in the next 
few months. The r decay particle distributions depend on the r spin direction. Thus the 
theoretical information on the r polarization will be valuable, since the expected number of 

events will not be large 0, lEI^- The study of r polarization is also needed, for instance, 
in i/^ — > Vf. oscillation experiments to disent ang le (z/g, e) events from background electron 
production following the — > oscillation [lO|. 

The study of lepton polarization in (i/^, V) and reactions is also of theoretical interest 
since the polarization observables may display peculiar sensitivities to the ingredients of the 
nuclear and reaction models, different to the ones shown by the cross sections. However 
the optimal neutrino-energy regime to obtain partially polarized leptons in non-longitudinal 
directions is limited. The reason is that for high energy (compared to it mass) the final 
leptons are 100% polarized with negative helicity. This is the case for the electrons in (z/g, e) 
reactions for most of the energies involved in experiments. In the case of muon production for 
some moderate energies, a non negligible transverse polarization component, though small, 
could be observed, and some examples will be shown below in this paper. More interesting 
is the aforementioned case of tau leptons due to the large value of its mass. This lepton will 
be the main focus of this paper. 

Previous studies of lepton polarization observables in neutrino induced reactions focus 
mainly on r polarization from nucleon targets for a range of kinematics of interest [lo| jisj]. 



2 



Neutrino detectors are based on neutrino-nucleus interactions (such as Argon in the ICARUS 
experiment), and a priori one would need to evaluate the lepton polarization observables 
in actual finite nuclei models. Genuine nuclear physics effects are usually neglected in 
these kind of reactions, since relatively high energies are involved. However for low to 
intermediate energy transfer nuclear effects do play a role in the inclusive cross section of 
(z/, I) reactions and in general cannot be neglected (Pauli blocking, finite size effects, long and 
short-range correlations, final-state interactions, sub-nuclear degrees of freedom, etc.). This 
is well known from electron scattering studies, that usually are the starting point of neutrino 
scattering models. For the same kinematics the (e, e') reaction, and the {v, I) reaction in the 
weak-vector sector are complementary to each other and make use of the same dynamical 
ingredients. Only the weak-axial current contribution makes a difference in the description 
of both reactions, since it can excite different nuclear modes. 

Only if a given model is able to reproduce to some extent the available (e, e') data, its 
predictions for neutrino reactions can be considered reliable. Different theoretical approaches 
along these lines on neutrino reactions in nuclei can be found among the recent literature 

available 0, H H El H H H H H 113 . 



In this paper we describe the {i^ul) reaction within the many-body framework of refs. 
|26[ I27I. Our model is an extension of previous studies on electron [28], photon jl^, and 
pion [30, 1^ interactions in nuclei. It is based on a Local Fermi Gas (LFG) description of 
the nucleus, accounting for Pauli blocking. The experimental Q-values are used to enforce 
a correct energy balance of the reaction. Additional nuclear effects, essential to describe e, 
7 and TT reactions, are built on top of the model, in particular: 

1. Coulomb distortion of the charged leptons, 

2. Medium polarization, which is taken into account through the random-phase approx- 
imation (RPA), 

3. Final-state interaction (FSI). 

The model was first applied to neutrino reactions in Ref. [2^, providing one of the best 
existing simultaneous description of inclusive muon capture, (z/^, yU~) and (z/g, e") reactions in 
^^C near threshold. Inclusive muon capture from other nuclei was also successfully described 
by the model. Apart from the description of the absorption of real photons by nuclei [29], 
the model describes rather well the (e, e') inclusive cross section of ^^C, ^''Ca and ^°^Pb for 
different kinematics, not only in the QE region, but also when extended to the A— peak and 
the dip region [i^. Recently the model has been extended to the neutral current sector, 
and one-nucleon knock-out reactions have been studied both for CC and neutral current 
driven processes j^]. A special effort has also been paid in Ref. |27[ to reliably estimate the 
theoretical uncertainty of our model. 

In this paper we extend the above model to the calculation of lepton polarization compo- 
nents in CC neutrino reactions. To our knowledge, only one previous study exists analyzing 
nuclear effects on lepton polarization [3^. This model basically consists on a relativistic 
local Fermi gas (LFG) including nuclear dynamical corrections such as some kind of rela- 
tivistic RPA correlations and an effective mass for the nucleons inside the nucleus (this is an 
approximate way of taking into account the FSI). However the model of [s^, first presented 
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3^ 3^, has not been tested in other nuclear reactions such as (e, e' 



In Sect. 2 we introduce the formalism and discuss the kinematics of the (z/;,/) reaction 
for muon and tau leptons in order to identify the most interesting cases. In Sect. 3 we 
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present results for differential cross sections and lepton polarization components, and draw 
our conclusions. 



II. FORMALISM AND KINEMATICS 

In this work we will study the inclusive neutrino-induced reaction depicted in Fig. 1. A 
neutrino Ui (or anti-neutrino Vj) with four momentum K^^ = {E^, k) exchanges a W boson 
with an atomic nucleus with initial momentum = (Mj,0), and a polarized lepton l~ 
(or /+) is detected with four-momentum K''^ = {E[^k'). In the inclusive reaction the final 
hadronic state is not detected. In this work we deal with the quasielastic channel, where 
the main contribution is due to one-nucleon emission. 




FIG. 1: Kinematics of the reaction, where we show the scattering plane and the directions of the 
two lepton polarization components P; and Pt- 



We write the unpolarized differential cross section as |2| 

^ _ £a^i \k'\G'^M, 



dVt'dE[ 27r2 



(1) 



where G is the Fermi weak coupling constant, f2' is the solid angle of the final lepton, and 
the quantity F has been defined as 

(Tfl^ \ — * -* 
2Wi + J^WA {E[ - \k'\ cosO) + W2{E[ + \k'\ cos9) 
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- W^5^ T ^ {e^EI + \kf - {E^ + El)\k'\ cose) (2) 

resulting from the usual contraction between the leptonic and hadronic tensors. Here 6 is the 
angle between k and k', and the =F sign in the last term correspond to the case of neutrino 
or anti-neutrino scattering. Finally the structure functions Wi are defined in the hadronic 
tensor as 
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FIG. 2: Allowed kinematics for (f^,//) reaction in the (a;,g) plane. The dotted line is the center of 
the quasielastic peak uj = {q^ — uj'^)/2mi^. The other lines show in each panel the value of g as a 
function of lo for a fixed value of the neutrino energy, and for several values of the scattering angle 0. 
Starting with the solid line for 6* = 0, the angles from down to up correspond to 6* = 10°, 20°, 30°, 40° 
and 50°. 



with eoi23 = +1 and the metric g^" = (+, — , — , — ). 

We assume that the final lepton polarization is measured in the direction defined by the 
vector s'^ verifying = —1. The polarized differential cross section (^J) can be written as 

S = (1 + s,Pn (4) 

where Eq, given in Eq. (1), is the cross section corresponding to unpolarized leptons and 
is the polarization vector. The relevant components of the polarization vector in the 
laboratory system are defined in Fig. 1 and denoted Pi (longitudinal, in the direction 
of A;'), and Pt (transverse to k' and contained in the scattering plane). Working within 
the standard model, as we do in the present work, it can be shown that the polarization 
component perpendicular to the scattering plane (Fig. 1) is zero. 0,11^. The expressions 
for the two polarization components. Pi and Pt, in terms of the hadronic structure functions 
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Wi are as follows 

[2Wi - j^W^j {\k'\ - E'lCosO) + W2{\k'\ + E^'cos^) - W^s^cos^ (5) 
T^((i^. + El)\k'\ - {E^E[ + \kf) cos^)| /F 

Pt = Tmi sin 9 (2W, -W2- + W,^ tWs^^ /F (6) 

Note that the transverse polarization Pt is proportional to the lepton mass. The larger values 
of this polarization component are then expected for tau leptons, while it will be negligible 
for electrons at the intermediate energies of interest for neutrino reactions. 





100 200 300 400 500 100 200 300 400 500 

uj [MeV] uj [MeV] 



FIG. 3: Allowed kinematics for {i't,t) reaction in the {LO,q) plane. The dotted line is the center of 
the quasielastic peak u; = (g^ — Lo'^)/2miy. The other lines show the value of g as a function of u 
for several fixed values of the scattering angle 6. Starting with the solid line for ^ = 0, the angles 
from down to up correspond to 6' = 2°, 4°, 6°, 8°, 10° and 12°. 



Let us begin with a discussion on the kinematics of the reaction in order to identify the 
cases of interest. For a fixed value of the neutrino energy, E^,, and from the definition of the 
four-momentum transfer, Q'^ = — K'^^ = {u},q), one easily obtains the modulus of the 
three-momentum transfer q = \q\, 

q'^^uj'^-mf + 2E^{E^ - u) - 2E^^ {E^ - u)^ - mf cos 9. (7) 

For a fixed value of the scattering angle, the above equation gives g as a function of cu, 
while the case ^ = leads to a parametrized curve corresponding to the the boundary of the 
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allowed kincmatical region in the (a;, q) plane. The identification of this region will help in the 
present work because we deal with a model describing the reaction around the quasielastic 
(QE) peak defined by a; = (g^ — u'^)/2m]\f, with m^v the nucleon mass. Therefore we must 
choose kinematics for which the QE peak lies inside the allowed region. Some examples are 
shown in Figs. 2 and 3 for muon and tau leptons, respectively. 

In Fig. 2 the muon case is displayed for four values of the neutrino energy ranging from 
El, = 500 to 1500 MeV. In each panel we show the possible kinematics in the {lu, g)-plane 
for several values of the scattering angle, ^ = 10°, 20°, 30°, 40° and 50°. The lower sohd lines 
correspond to the limit value ^ = 0, setting the lower boundary of the allowed kinematical 
region. The maximum value of the energy transfer uj^na^ = Ey — mi sets the right end of 
the boundary. In the same plots we also show the position of the maximum of the QE peak 
uj = {(f — u}'^)/2m]si with dotted lines. Note that in all cases the QE peak region is inside 
the allowed region for {v^^ij) reactions. In the plots we only show the low to intermediate 
energy transfer part uj < 500 MeV since we are interested in the non relativistic regime 
where our nuclear model is safely applicable. This also means that the momentum transfer 
q must be below 600 or 700 MeV/c, since relativistic corrections at these values start to be 
important. Although our nuclear model is based on a fully relativistic description of the 
Fermi gas, the RPA and FSI corrections are non relativistic. Therefore for this work we are 
forced to kinematics with low to intermediate energy and momentum transfer. The plots in 
Figs 2 and 3 arc very useful to this end since they give us a clear picture of the kinematical 
changes when the scattering angle is increased. 

The case of r leptons is considered in Fig. 3 for neutrino energies going from E^ — A 
to 15 GcV. This time the scattering angles shown are 6 = 0, 2°, 4°, 6°, 8°, 10° and 12°. In 
contrast to the muon case, for tau leptons the QE region is forbidden for some kinematics. 
In particular for low neutrino energy and low scattering angles. For instance, for Ei,=4: GeV 
the maximum of QE peak is below the ^ = curve for uj smaller than 200 MeV. In order 
to cover this region, one should go to momentum transfers well above 800 MeV/c, where 
relativity would play a role. Therefore we must go to larger neutrino energies, above ~ 7 
GcV, in order to cover the QE region for small values of q and uj. Moreover, for large values 
of Ey the momentum transfer strongly increases with the scattering angle, and then we force 
it to fall in the few degrees region to guarantee non relativistic kinematics. 

III. RESULTS AND CONCLUDING REMARKS. 

Our results showing the effects of the different corrections implemented into a nuclear 
model are summarized in Figs. 4-8. With dotted lines we show the results of our model 
without RPA and without FSI corrections. Hence these results account basically for the 
impulse approximation with neutrino-nucleon interaction oiV — A type. We use the Galster 
parameterization for the nucleon vector form factors and a dipole dependence for the axial 
form factor, while the pseudo-scalar form factor is related to the later by the partially 
conserved axial current (PCAC) hypothesis. These results include also the following effects: 

1. Pauli blocking, through a LFG description of the nucleus. This implies an additional 
dependence on the experimental nuclear density. 

2. Correct energy balance of the reaction using experimental Q— values. In the figures 
we show g as a function of lu minus the Q value. 
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FIG. 4: Differential cross section and polarization components for the Ar(f^,;u) reaction as 
a function of the energy transfer minus the experimental Q-value, and for three values of the 
momentum transfer q = 300, 400, and 500 MeV/c. The incident neutrino energy is fixed to 500 
MeV. 



3. Coulomb distortion of the charged leptons. 

Ivlore details on the parameters of CC current and on the model are given in our previous 
works [26, 2;^. 

In the same plots we show with dashed lines the results including RPA corrections, that 
take into account modification of the nuclear medium through interactions between particle- 
hole and A— hole excitations. We use an effective nucleon-nucleon interaction, with pion and 
rho exchange in the vector-isovector channel, and corrections due to short-range correlations. 
We refer the reader to Ref. for details about the RPA series and for explicit expressions 
of the nuclear structure functions Wi entering in the CC neutrino cross section. 

Finally we show with solid lines the results with the full model, including also the FSI. We 
account for relevant reaction mechanisms where two nucleons participate, by dressing the 
nucleon propagators in the nuclear medium. In particular these effects change the dispersion 
relation of the nucleon, which in some works is taken into account by the over-simplified 
method of introducing an effective mass for the nucleon [s^ . 

We start the discussion with Fig. 4 where we show typical results for the (z/^, fi) reaction. 
In the figure we show the relevant observables, namely the differential cross section and the 
two polarization components. Pi and Pt in the laboratory system. The neutrino energy is 
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FIG. 5: The same as Fig. 0]for tau neutrinos of 7 GeV. 



500 MeV, and three values of the momentum transfer q = 300, 400 and 500 MeV/c are 
considered. These values of the momentum transfer correspond to lepton scattering angles 
of 30° and above. The cross section shows the typical behavior of the QE peak and we 
can see how the different nuclear effects modify this observable. First the RPA produces 
a big reduction and small shift of the peak (compare the dotted to the dashed lines). The 
quenching of the cross section is large for low momentum transfer (almost a 50% reduction 
at the peak for q = 300 MeV/c) and diminishes with q (below 15% for q = 500 MeV/c). On 
the other hand the FSI produces a further reduction of the RPA results and an important 
enhancement for high energy transfer. This a consequence of the modification of the nucleon 
dispersion relation in the medium through the dressed nucleon propagator inside the nucleus. 

Concerning the polarization observables, we see that the longitudinal component Pi is 
very close to —1 for the three kinematics, while there is a small but appreciable transverse 
component (around —0.2), quite independent of q. The inclusion of the RPA does not 
change these results. The reason is that the polarization components are obtained as a ratio 
between linear combinations of nuclear structure functions and the RPA changes similarly 
numerator and denominator. The same can be said for the FSI effects, except for the high-u; 
region. The LFG is unable to contribute to the high energy tail of the cross section, while 
the model with FSI is able to describe this region where, however, the QE cross section is 
rather small. 
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FIG. 6: The same as Fig. El for tau neutrinos of 10 GeV. 



We have generated results for higher values of the muon-neutrino energy corresponding 
to the kinematics of Fig. 2. However in all the remaining cases the transverse polarization 
component computed is negligible and the muon can be considered polarized with negative 
helicity. 

An example of the results found for tau leptons is displayed in Fig. 5, where we show 
the {vri t) cross section and polarization observables for the same values of q as in Fig. 4, 
but this time for = 1 GeV. By inspection of Fig. 3 we see that for the three values of 
q considered the maximum of the QE lies inside the allowed kinematical region. The case 
q = 300 MeV/c is closer to the boundary, since the cross section suddenly ends soon above 
the maximum, while for higher values of the momentum transfer the allowed region extends 
to higher energy transfer. The effects seen here over the cross section due to RPA and FSI 
are similar to the muon case studied above, with the exception of the missing high energy 
tail that lies now in the forbidden region. 

Concerning the polarization components, the longitudinal one is still negative but now 
is well above —1. On the other hand, an appreciable Pt component appears, taking values 
ranging from ~ —0.2 to ~ —0.25. Thus the magnitude of Pt increases with q, and, as in the 
muon case, the RPA and FSI effects are rather small on the polarization observables. 

A more clear scene with similar results can be seen in Fig. 6 for Ei, = 10 GeV, corre- 
sponding to the kinematics of the third panel in Fig. 3. In this case the QE peak is well 
inside the allowed region for the three values of q. The polarization component Pi is now 
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FIG. 7: The same as Fig. [3 fixing the scattering angle 6 instead of the momentum transfer q. 
Three values of ^ = 0°, 2° and 4° are shown. 



closer to —1, while the magnitude of Pt is smaller than in the former case. Results for 
El, = 15 GeV (not shown) indicate that Pi almost reaches the limit value ~ —1, while Pi 
becomes very small for high neutrino energy, as expected. 

Another kind of plot of interest for ongoing neutrino experiments is shown in Fig. 7 for r 
neutrinos with energy of 7 GeV. This time we fix the scattering angle instead of the momen- 
tum transfer. Therefore by changing u we are always inside the allowed kinematical region, 
running along some of the curves shown in the second panel of Fig. 3, never crossing the 
boundary, and the cross section is defined for every value of u shown in the plots. Results 
are shown for small angles, 6 = 0, 2°, and 4°, in order to reach not too high values of the mo- 
mentum transfer. For these angles the maximum of the QE peak is crossed at g ~ 300, 400, 
and 600, respectively. Except for ^ = 4° at the high u tail where the g-values are perhaps 
too high, and a word of caution is needed since some important relativistic corrections are 
expected, we can safely trust the results for lower cu- values where q is moderately low. Again 
we can see the important reduction of the cross section due to RPA correlations for low u. 
The magnitude of this reduction decreases with the scattering angle. The FSI produces a 
reduction for low u and an increase for high values of the energy transfer. The polarization 
component Pi is more or less independent on the angle and takes values between -0.9 and 
-0.8. The transverse polarization Pt is zero for = by definition, Eq. (6), while it is not 
negligible at all for the remaining angles, taking values between ~ —0.2 and -0.5. Concern- 
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FIG. 8: Modulus and angle of the polarization vector with respect to the — / direction for tau 
neutrinos of energy 7 GeV, and for three scattering angles, corresponding to the kinematics of fig. 

m 



ing the RPA and FSI corrections on these polarization observables, the effect is found to be 
neghgible for low energy. However at the high energy tail and for 6 = the net RPA+FSI 
effect makes Pi to change from ~ —0.6 to ~ —0.8. A less important, but still appreciable 
change is found in Pt for 9 = 2°. For higher values of the scattering angle the effect is again 
neghgible. 

To end the discussion, in Fig. 8 we show another representation of polarization observ- 
ables, namely, the total polarization, defined as the modulus of the polarization vector 

1^1 = ^|pf+Pl 

and the angle with respect to minus the longitudinal direction {—k'), 

e = arctan(Pt/P;). 

— > 

The total polarization \P\ takes values between 0.8 and 0.9 for 9 = and increases with 
the scattering angle. For 9 = the angle B is also zero, meaning that the polarization 
vector points to the —k' direction, without transverse components. In this case the total 
polarization has the meaning of fraction of particles with negative helicity. A value different 
from unity means that the interaction with the nuclear target produces a small fraction of 
tau leptons with positive helicity. A value close to one means that a large percentage of the 
leptons exit with the spin pointing to the same direction. The most probable direction of 
the lepton spin is determined by the angle G. For 9 = 2° this angle is between ~ 15 and 25°, 
and slightly increase for 9 = 4°. The results of Fig. 8 are showing a non negligible increase 
of the total polarization due to RPA+FSI effects, for 6' = 0. These effects are negligible 
over this observable for 9 = 2 and 4°. On the other hand, for 9 = 2° we find an appreciable 
reduction of the polarization angle 6 due to RPA+FSI, that again is negligible for higher 
scattering angles. 
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Summarizing this work, we have computed the cross section and polarization observables 
for neutrino induced CC reactions in nuclei at the QE peak. We have focused this work 
on the case of r leptons of interest for neutrino oscillation experiments. The RPA plus FSI 
nuclear effects have been evaluated for intermediate energy and momentum transfer. These 
corrections are essential in the cross section, especially for low energy, but are partially 
reduced, due to cancellations, on the polarization observables. However, we have identified 
some particular kinematics, at very low scattering angles, where these nuclear corrections 
are of some importance to determine the correct magnitude and angle of the polarization 
vector. 
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